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We show that a general two-Higgs-doublet model (THDM) with a very light CP-odd scalar 
(A) can be compatible with p parameter, Br(6 — > 37), Rt, Ah, (g — 2) M of muon, Br(T — > A-y), 
and the direct search via the Yukawa process at LEP. For its mass around 0.2 GeV, the muon 
(g — 2) M and Br(T — > Af) data require tan/3 to be about 1. Consequently, A can behave like a 
fcrmiophobic CP-odd scalar and predominantly decay into a 77 pair, which registers in detectors of 
high energy collider experiments as a single photon signature when the momentum of A is large. We 
compute the partial decay width of Z — > AAA and the production rate of // — > ZAA — » Z + "77", 
/'/ -> W ± AA -> W ± + "77" and // -> H+H~ -> W + W~AA -> W+W~ + "77" at high energy 
colliders such as LEP, Tevatron, LHC, and future Linear Colliders. Other production mechanisms 
of a light A, such as gg — > h — ► AA — > "77", are also discussed. 

PACS numbers: 12.60.Fr, 14.80.Cp, 12.15. Ji 



I. INTRODUCTION 



One of the main tasks of the current and future high energy colliders is to find the Higgs boson of the Standard 
Model (SM), or some other scalar particle(s), if there is any, predicted by the extensions of the SM. The mass spectrum 
of these scalars as well as their decay channels depend on the assumed model. Recently, the possibility of a Higgs 
boson decaying into a pair of light CP-odd scalars was considered in Ref. Q. As pointed out in that paper, a very 
light CP-odd scalar (A) can arise in some extensions of the SM, such as the minimal composite Higgs model or 
the next-to-minimal supersymmetric model Q. An interesting aspect of a light A particle is that if its mass (Ma) 
is less than twice of the muon mass (m^), i.e. less than about 0.2 GeV, it can only decay into a pair of electrons 
(A — > e + e~) or photons (A — + 77). Hence, the decay branching ratio Br(A — ► 77) can be sizable. Furthermore, at 
high energy colliders, the light CP-odd scalar A can be produced with such a large velocity that the two photons 
from its decay are highly boosted and seen by the detector as one single-photon signature. Therefore, the production 
of a pair of AA is identified by the detector as a pair of photons when each A decays into its diphoton mode. The 
subsequent signature as a diphoton resonance (e.g., in the case that a Higgs boson decays into a pair of A) or photon 
cascades (e.g., in the case that a A particle is radiated from a fermion line) provides an interesting window for the 
experimental search of the scalar particles that may be responsible for the breaking of the electroweak symmetry . In 
general, we expect the decay branching ratio of A — > 77 to decrease rapidly when the di-muon channel (A — ► /i + /i~) 
becomes available as Ma increases. This is because the Yukawa coupling of A-/.i + -/j,~ is larger than that of A-e + -e~ 
by the mass ratio (m M /m e ) ~ 200. Since we are interested in the phenomenology of having a light A decaying into 
a pair of photons, we will restrict our discussion for Ma less than about 0.2 GeV, though in principle, that is not 
necessary as long as the decay branching ratio Br(^4 — > 77) is not too small to be observed experimentally. 

A simple extension of the SM is the two-Higgs-doublet model (THDM) ||, which has been extensively studied 
theoretically and experimentally. For example, a constraint on the mass of the charged Higgs boson H + in the 
THDM was carefully examined at the CERN LEP as a function of its decay branching ratios into the t + v t and cs 
modes B . Studies on searching for a light A in its associated production with a bottom quark pair were also done M , 
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as a function of Ma and tan (3 (the ratio of the two vacuum expectation values of the Higgs doublets in the THDM) . 
Though some useful constraints have been obtained by the LEP experiments, we show in this work that a very light 
A (with Ma < 2m M ) is still allowed in the THDM. This low Ma value, in the context of a THDM, can induce large 
contributions to the p-parameter unless other parameters of the model adjust to counteract this effect. In particular, 
as to be discussed later, the masses of the charged scalar (H + ) and the heavy CP-even Higgs boson (H) have to be 
approximately equal. Also, the value of the mixing angle a has to be such that cos(/3 — a) becomes small in order to 
suppress the Higgs boson contribution to the p-parameter. 

Another important experimental data to constrain a light A in the THDM is the measurement of the muon anoma- 
lous magnetic moment. The recent measurement at BNL || strongly disfavors such a model when compared with 
certain theory calculation Jt],||]. Nevertheless, other theory calculations of the SM contributions to the muon magnetic 
moment show a better agreement with the BNL data || . Consequently, a light A in the THDM can still be compatible 
with data though the parameter space of such a THDM is tightly constrained. 

In the next section we will examine all the relevant low energy data (including the (g — 2) M of muon, Br(T — > A7), 
p-parameter, b — ► 57, Rb and Ab) to determine the allowed parameter space of the THDM with a very light A. In 
section III, we consider the decay widths and the decay branching ratios of every Higgs boson predicted in this model. 
In section IV, we study the potential of the light A boson as a source of the distinctive photon signal at colliders. It 
can happen in either the decay mode of the neutral gauge boson Z — > AAA — > "777", or the production processes of 
Higgs bosons, such as // -> ZAA -> Z + "77", /'/ -> W ± AA -> W ± + "77" and // -► H+H~ -> W+W~AA -> 
W + W~ + "77". Section V contains our conclusion. 



II. CONSTRAINTS FROM LOW ENERGY DATA 



In the THDM, the couplings of fermions to Higgs bosons are proportional to the fermion masses. In the type-I of the 
THDM, only one of the two Higgs doublets couple to fermions via Yukawa couplings, and in the type-II of the THDM, 
one of the Higgs doublets couples to the up-type fermions (with weak isospin equal to 1/2) and another couples to the 
down-type fermions (with weak isospin equal to —1/2). Hence, the couplings of Higgs bosons to fermions generally 
depend on the value of tan (3. In case that the coupling of the Higgs bosons to fermions is large and the mass of the 
Higgs boson is small, the radiative correction to the low energy data can be sensitive to the Yukawa interactions. 
Hence, the low energy data can be used to impose important constraints on the masses and the couplings of the Higgs 
bosons. To examine the allowed range of Ma and tan/3 in the THDM, we shall consider the precision data on the 
anomalous magnetic moment of the muon, the decay branching ratio of T — > Aj, b — > S7, the p-parameter, the decay 
branching ratio of Z — > bb, i.e. Rb, and the bottom quark asymmetry Ab measured at the Z-pole. 



A. Constraint on tan/3 from the (g — 2) M of muon 

The magnetic moment of muon is defined as 

eh 

flf, = (1 + cfy) 



2m M 



with a, = , (1) 

where a M is the muon anomalous magnetic moment, which is induced from radiative corrections. The SM prediction 
includes the QED, weak and hadronic contributions. Among them, the hadronic contribution has the largest uncer- 
tainty, and the bulk of the theoretical error is dominated by the hadronic vacuum polarization (h.v.p.). There are 
a number of evaluations of the h.v.p. corrections; four recent results were extensively discussed in Ref. Q. After 
comparing various theory model predictions with the precise experimental data, which is 

10 11 x (experiment) = 116 592 023 ± 140 ± 60 , 

Yndurain concluded that the discrepancies between the world averaged experiment data (exp.) and the theory pre- 
diction of the SM contribution (theo.) are jj): 

10 11 x Aa M (exp. - theo.) = 422 ± 152(exp.) ± 77(theo.) (DH) 
10 11 x Aa M (exp. - theo.) = 358 ± 152(exp.) ± 112(theo.) (J) 
10 11 x Aa M (exp - theo.) = 233 ± 152(exp.) ± 104(theo.) (AY) 
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10 11 x Aa M (exp. - theo.) = 119 ± 152(exp.) ± 115(theo.) (CLY) (2) 

In the above result, DH stands for the analysis of Ref. J for that, as yet unpublished, of F. Jegerlehner, quoted 
in Ref. 0, AY indicates the result of Ref. pi] , and CLY is the 'old' result of Ref. jl2| after being corrected for the 
new favored value of higher order hadronic corrections |J . 

In the general THDM, as well as in supersymmetric models, et M can receive radiative corrections at the one loop 
level from the couplings of A, h and H to muons in triangle diagrams (see Appendix |a|) J|,|l3],[l4]] . As expected, 
the size of the radiative corrections is proportional to the coupling of muon to Higgs bosons. Moreover, the loop 
integral reaches its maximal value when the mass of the scalar boson in the loop becomes negligible, and diminishes 
as the scalar mass increases. It was concluded in Ref. || that a light CP-even Higgs boson (h) can be responsible for 
the apparent deviation of the BNL measurement of a M from the DH prediction of the SM contribution, at the 90% 
confidence level ( CL), in the framework of a type-II THDM, in which the other Higgs bosons are heavy (of the order 
of 100 GeV). It was found that the model parameters have to satisfy the following requirements: 

m T ^ M h < 2m B , 
sin(/3 — a) ~ , 

30<tan/3<35. (3) 

In the case that only A is light and the other Higgs bosons are heavy, the one-loop contribution to a M is negative.^ 
Therefore, this type of new physics is strongly disfavored according to the theory prediction of the SM contribution 
provided by DH; however, it can still be compatible with the other SM theory calculations (cf. Eq. ||). Furthermore, 
it was found in Ref. fl5|| that a two- loop contribution to can be sizable when A is light. As compared to the 
one-loop graph, a two-loop graph can contain a heavy fermion loop. The Yukawa coupling of the heavy fermion 
(with mass m/) in the second loop together with the mass insertion of the heavy fermion will give rise to (m//m M ) 2 
enhancement which can overcome the extra loop suppression factor of 1/167T 2 . Because the two- loop contribution can 
be even larger than the one-loop contribution, the contribution from a light A to can become positive when Ma is 
not too small. Hence, in general, a two-loop calculation of a light A contribution to the muon magnetic moment a M 
yields a better agreement with the experimental data than a one-loop calculation. For that reason, in the following 
numerical analysis, we shall apply the two-loop calculation presented in Ref. jl5| to test the compatibility of a light 
A THDM to the experimental data (cf. Eq. |^). For completeness, we summarize the relevant formula in Appendix ^ 
to clarify the contributions included in our numerical analysis. 

Although we are using a two-loop calculation for our numerical analysis, it is useful to examine a few features 
predicted by the one-loop calculation. Since a one-loop contribution to from a light A is always negative and the 
central value of a Al (exp. - theo.) is positive, the potentially large loop contribution to a p has to be suppressed by a 
small Yukawa coupling in order for the model to be compatible with data. For a type-II THDM, this implies a very 
stringent bound on tan/3, because the coupling of muon to A is directly proportional to tan/3. (The coupling strength 
of A-fj, + -fi" in the type-II model is (m M /i>) tan/3, where v is the weak scale, ~ 246 GeV.) In particular, assuming the 
CLY prediction for the SM contribution and applying the two-loop calculation to include the light A contribution, we 
find that there exists an upper bound on tan/3 at the 95% CL. For a 0.2 GeV pseudo-scalar, 

tan/3<2.6. (4) 

This new bound is stronger by a factor of two than a previous one jl6| , which was obtained from the one loop 
contribution, together with the old experimental data with an error of 84 x 10 -10 in the measurement of (g — 2) M . 
In Fig. [l], we show the regions in the tan/3 versus Ma plane allowed by the data at the 95% CL. Three different 
curves are displayed depending whether the SM prediction is given by the CLY, AY or J calculation. There is no 
allowed region according to the DH calculation. For the type-II THDM, the allowed regions are below the curves, for 
the type-I THDM, they are above the curves.^ For completeness, we also show in Figs. ^ and || the allowed regions 
in the type-I and type-II THDM, respectively, for a wide range of Ma and tan/3. In the above figures we did not 
show the constraints derived from the CUSB Collaboration search for T — > Aj at CESR |Q. For that, we refer the 
readers to Ref. [|J which has an extensive discussion on this data to constrain a light CP-even or CP-odd scalar in 
the THDM. As noted there, various theory analyses indicated that the high order QCD corrections to this decay rate 



1 In contrast, the one- loop contribution of a light h to a M is positive. 

2 At the one- loop level, the pseudoscalar contribution to (g — 2) M in a type-I model can be obtained from that in a type-II 
model after replacing tan j3 by cot 0. 
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can be large. Because of that, we shall consider in this paper tan/3 around 1 to be consistent with the CUSB data 
for a 0.2 GeV pseudoscalar. Specifically, we shall take tan/3 to be either 0.5 or 2 in our following discussions. 

Finally, we note that a similar constraint on tanj? can be obtained from examining the production of e~e + — > bbA 
at LEP energies ||. However, as shown in Ref. Jl6|], the (<? — 2) M data gives a more stringent constraint than that 
obtained from direct searches of Higgs bosons via the above Yukawa process. 

B. Constraint on M H + from the decay of b — *• S7 

The previous (g — 2) M analysis is only sensitive to a light pseudoscalar A when the masses of the CP-even scalars h, 
H and the charged Higgs boson are all of the order of hundred GeV. To further constrain the parameter space of 
the THDM, we now turn our attention to a low energy observable that is sensitive to a charged Higgs boson in case 
that tan j3 is not large. That is the rare decay process b — > 37. 

For the hadronic flavor changing neutral current decay process b — > sj, new physics effects at the weak scale can 
be parameterized by the couplings (Wilson coefficients) of an effective Hamiltonian |ll| . The only scalar contribution 
to the Wilson coefficient C7 in a THDM comes from the standard penguin diagram, where the charged Higgs scalar 
H ± couples to top and bottom, and then to top and strange quarks. 

The H ± tb coupling in the type-II THDM is given by 

C = V2V tb ( — tan t L b R + — cot f3 t R b L ) H + + h.c. (5) 

where Vtb is the Cabbibo-Kobayashi-Maskawa (CKM) matrix element. The .ff is coupling is defined similarly with 
the appropriate substitution b — > s. In the type-I model, the factor tan/3 is substituted by — cot/3. For a type-II 
model, a value of tan/3 < 1 will induce a large i,Rbi J H + coupling, which will strongly modify the Wilson coefficient 
C7 and then increase the predicted b — > sj rate. This can be alleviated only if the charged H + is massive enough. 

Constraints on tan/3 versus Mjj+ have been obtained by F. Borzumati and C. Greub in the first reference of [H| 
for a few possible experimental upper bounds on Br(_B — * X s j), ranging from 3.0 x 10~ 4 to 4.5 x 10~ 4 . Currently, 
the reported experimental measurements are 

Br{B X sl ) x I0 4 = 3.15 ± 0.35 stat ± 0.32 sys ± 0.26 mode , by CLEO [20], (6) 
Br(B -> X sl ) x I0 4 = 3.II ± 0.80 stat ± 0.72 sys by ALEPH fl|, (7) 
Br(B-X s7 ) x I0 4 = 3.36 ±0.53 stQt ± 0.42^54) ™^ by BELLE |§. (8) 

In this study, we will quote the result given in Ref. pi| for an upper bound of 4.5xl0~ 4 on Br(£? — > X s j). It was 
found that, for example, if tan/3 = 0.4, the mass of must be larger than about 300 GeV. On the other hand, 
we will not use the detailed information on the lower bounds of M#+ for tan/3 < 1, because in that case we would 
rather use the stronger bounds obtained from examining the bb decay rate of the Z-boson, i.e. (see next section). 
When tan/3 is much greater than 1, the small effect from the tnbLH + coupling is compensated by the large effect 
from the ti,b]iH + coupling; in such a way that the Br(fe — > sj) prediction tends to stay at a certain minimal value. 
This minimum only requires Mjj+ > 165 GeV for tan/3 3> I. In summary, for tan/3 > 1, the Br(& — > sj) data requires 
M H + > 165 to 200 GeV at the 95% CL. 

The situation for type-I models is different. In this case, both the tnbLH + and tLbnH + couplings decrease for 
values of tan/3 greater than 1. This does not mean that M H + can become arbitrarily small. There is an unstable 
behavior (due to the large scale dependence) for the prediction of Br(6 — > sj) when M H + is less than about 100 GeV 
and tan/3 ~ 1 |l9|. This problem of unstability gets worse when tan/3 < 1. 

C. Constraint on Mh — M H + from Ap 

The effect of scalar fields (A, h, H and H + ) on Ap have been reported for a general THDM in the literature [p3[-p5|. 
When Ma becomes negligible, the contribution to Ap can grow quadratically with the masses of the other scalars 
(See the appendix of Ref. |2q|). The only way to keep Ap small is to have cancellations among potentially large loop 
contributions from each Higgs boson. For example, as sin 2 (/3 — a) — 1 , this cancellation can take place between the 
contributions from H ± and H, and requires a certain correlation between their masses, (a is the mixing parameter 
between the two CP-even Higgs bosons h and H.) As shown in the appendix of Ref. j25|, this correlation depends on 
the value of the coefficient sin ifi-a). In Fig. I we show the different allowed regions in the M#+ versus Mh plane 
for three different values of sin 2 (/3 — a). These regions or bands correspond to the 95 % CL limits for Ap |25| with 
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- 1.7 < Ap x 10 3 < 2.7. 



(9) 



For simplicity, we shall assume Mh ~ M H + and sin 2 (/3 — a) = 1 for our detailed numerical analysis. It should 
be noted that the allowed band for sin 2 (/3 — a) = 1 only depends on Mh and M H +, whereas the allowed bands for 
sin 2 (/3 — a) < 1 do depend on M/j. In Fig. || we have used a low value of Mh = 110 GeV. If we use a higher value, for 
instance Mh = 130 GeV, the bands keep the same width but slightly shift downwards (to a somewhat smaller slope). 

We have checked that within the allowed parameter space constrained by the p-parameter, the new physics con- 
tribution from the THDM with a light A to the S-parameter is typically small as compared to the current data: 
AS = — 0.07± 0.11 [p6| . Using the analytical result in Ref. p7f , we find that AS approximately equals to —0.02 when 
Mh and Mh+ are about the same, and reach —0.05 when Mh and Mh+ differ by a few hundred GeV. Although 
AS* can reach —0.1 when Mh and M H + differ by about a TeV, this choice of parameters is already excluded by the 
p-parameter measurement. Hence, we conclude that the current .^-parameter measurement does not further constrain 
this model. Needless to say that the above conclusion holds for both the type-I and type-II THDM at the one-loop 
order. 



D. Constraint on M H + versus tan/3 from Rb and At 

The contributions to the Z — > bb hadronic decay branching ratio (Rb) and the forward-backward asymmetry (Ab) 
of the bottom quark in Z decays are given in terms of the effective Zbb couplings p5| : 

C = — (g L b Ll n L + g R burfbn) Z^ , (10) 



1 l 

2 + r 



1 2 , x 

9R = + °9 R ■ 

Here, <5g's contain the SM as well as the new physics (THDM) contributions at one loop order: 

5g L = Sgf M + 5gT w = -0.4208 + 6g n L ew , (11) 
Sgn = Sg s B M + 6gl ew = 0.0774 + 6g% w , (12) 

where the SM values are for m t — 174 GeV and Mh = 100 GeV [] |^| . Because of the left-handed nature of the weak 
interaction, the value of gL is higher (by about 5 times) than g^. 

As shown previously, the (g — 2) p data requires tan/3 to be less than about 2.6 for a 0.2 GeV M^.Q In that case, 
the contribution from the neutral Higgs boson loops to Rb is negligible and the dominant contribution comes from the 
charged Higgs boson loop p5|| . For a small tan 0, the H + loop contribution comes mostly from the tpbLH + coupling, 
which is proportional to (mt/v) cot/3, cf. Eq. (|l3|). Since the type-I and type-II models coincide on this coupling, 
the bounds from the Rb measurement applies to both of the models. The H + loop contributions to 5gL and Sgn are 
given by 

Sg2T = cot/3) 2 f-g- - m °Z R \ , (13) 

yL 16tt 2 V v A J \R- 1 (R- 1) 2 y ' V ' 

X new 1 f m b . Q \ 2 ( R R\og R \ 

8gl ?«" = -__ — tan/3 - — --7= — , (14) 



16tt 2 Id W V-R- 1 !) 2 / 
with i? = mf/Mf I+ . Expanding Rb and to first order in Sg R ew and 5g ne ' w , we obtain 



= - 0.7785^£ etu + 0.1409(5^, (15) 
A 6 = Af M - 0.2984<5(7£ eto - 1.6234(5^ , (16) 



3 In the heavy top quark mass expansion, Rb depends on Mh through log(Mh/mz); therefore, the dependence of Rb on Mh is 
small. 

4 Unless specified otherwise, we shall assume the theory prediction by CLY in the following discussion. 
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with R^ M — 0.2158 and Af M — 0.935. As indicated in the above equation, At depends more on g^ than on g^; the 
opposite is true for Rb- For this reason, and because the experimental uncertainty of Rb is significantly smaller than 
that of Ab, the asymmetry At, is not nearly as effective to constrain the parameter region of the THDM than Rb when 
tan/3 is of the order of 1. Hence, in this work, we will only consider the constraints imposed by the Rb measurement, 
for which the experimental limits are p9|: 



The allowed region, at the 95% CL, in the tan/3 versus M H + plane is shown in Fig. [|, in which the constraints imposed 
by the (g — 2) p and the b — > S7 data are also included. It is interesting to compare our Rb bound (the lower solid 
curve) with the Fig. 6.1 of Ref. (same as the Fig. 5 of Ref. |24| ). There, the experimental central value used for 
Rb was 0.21680, with a smaller (lcr) error of ±0.00073. Since the new R^ xp is only 0.81<r above the SM value that 
we use here, the bounds on M H + become less stringent^]. 

In summary, Fig. || shows the allowed region in the tan (3 versus Mjj+ plane for the type-I and type-II models with 
a light (0.2 GeV) CP-odd scalar A. The (g — 2) M data imposes an upper bound on tan j3 for a type-II THDM, and the 
value of this upper bound depends on the model of theory calculations. Note that there is no DH curve in the figure. 
This is because a THDM with a very light A cannot be compatible with data according to the DH calculation of the 
SM contribution to a M . Similarly, the lower bound on tan /3 imposed by the (g — 2)„ data for a type-I THDM can be 
easily obtained from Fig. [fl. The lower bound from Rb holds for either type-I or type-II THDM, and is not sensitive to 
the actual value of the mass of the light A, because for tan/3 ~ 1 the charged Higgs boson loop dominates. A similar 
conclusion also holds for the constraint imposed by the b — > 57 data. Furthermore, the b — > sj data does not provide 
a useful constraint for a type-I model when tan/3 > 1. For clarity, we summarize the constraints from the relevant 
low energy data in Table | for each type of THDM. Although in the type-II THDM, the value of tan /3 is bounded 
from above to be less than about 2.6, it can take on any value above the lower bound imposed by the Rb data. On 
the other hand, the value of tan /3 in the type-I model cannot be too large because for a very large value of tan /3, the 
decay width of the lighter CP-even Higgs boson (h) can become as large as its mass, so that the model ceases to be 
a valid effective theory. 

In the next section, we shall discuss the decay and production of the Higgs bosons in the THDM with a light A. 
For simplicity, we shall only discuss the event rates predicted by a type-II model, and the value of tan /3 is taken to 
be either 0.5 or 2 to be consistent with all the low energy data discussed in this section. As shown in Fig. [|, in a 
type-II THDM with a light A, the masses of the Higgs bosons other than the light CP-even Higgs boson h have to be 
around 1 TeV when tan/3 = 0.5. Hence, its phenomenology can be very different from the model with tan/3 = 2, in 
which a few hundred GeV heavy Higgs boson H and charged Higgs boson H ± are allowed. 



In this section, we shall examine the decay branching ratios of the Higgs bosons predicted by the type-II THDMs 
with a very light CP-odd scalar A. To be consistent with the p-parameter analysis, the mixing angle a as well as the 
Higgs boson masses M H + and Mh have to be correlatively constrained, cf. Fig. ||. Without losing generality, in this 
section we shall assume sin(/3 — a) = 1 to simplify our discussion. 

As noted in the previous sections, a CP-odd scalar A, with its mass around 0.2 GeV, can only decay into a e~e + 
pair or a photon pair. Though the decay process A — > 77 can only occur at loop level, its partial decay width may 
compete with the tree level process A — > e~e + . This is because the mass of the electron is very tiny as compared to 
the electroweak scale v, and the partial decay width of A — > e~e + is suppressed by (m e /v) 2 . 

To clarify our point, we note that for a — f3 — ir/2, the couplings of Higgs bosons and fermions in the THDM are 
given by 



R» x p = 0.21648 ± 0.00075 . 



(17) 



III. DECAY BRANCHING RATIOS OF HIGGS BOSONS 



C lfq = -ffh - X ddH + Y uuH 

1) V V 

+ iX^d^dA + iY^-u^uA, 



(18) 



v v 



5 For our Rb analysis, we used the same input parameters as those given in Appendix H of Ref. J25|, but with an updated 
value for Rf xp . 
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where X = tan (3,Y = cot (3 for the type-II model, and X — — cot (3, Y — cot (3 for the type-I model. (Here, / stands 
for any fermion, u for an up-type fermion, and d for an down-type fermion.) Therefore, the partial decay widths of 
Higgs bosons into fermion pairs are: 



Nr. ~o , , / . . m2 f 



3/2 



/ 2 \ !/ 2 

r -//=£^ M 4 1_4 ^ij ' (20) 

where — ^-x(l, X, or Y ) is the coupling defined in £jj S , and iV c is the color factor (which is 3 for quarks 

and 1 for leptons). 

The partial decay width of A — > at the one-loop order arises from fermion loop contributions, which yield 

H^NtCftCfTfFiT,) (22) 
/ 

where r/ = (2m,f/MA) 2 , N[ — 1 (3) for leptons (quarks), Q/ and m/ are the electric charge (in units of e) and the 
mass of fermion, respectively. Also, Cf — cot (3 (tan /3) for up-type quarks (charged leptons and down-type quarks) in 
a type-II THDM. Furthermore, F(x) is given by 

■ ( arcsin -4= ) for x > 1 

F(x)=<| \ ^) (23) 

( arccosh^ — ] for x < 1, 



Wx 2 



In Figs. |^ and ^, we show the life-time (multiplied by the speed of light) of A and its partial decay branching ratio 
Br(A — * 77), as a function of tan/3 for various Ma values. As indicated, the typical life-time of a light CP-odd scalar 
A (with mass around 0.2 GeV and tan/3 around 1) is about meter, so the decay length of a 50 GeV A boson 

is about 0.25 meter. This unique feature of a light scalar A boson can be used to improve identifying such an event 
experimentally. For tan/3 ~ 1, about half of the time, the light A can decay into a photon pair, other than a e~e + 
pair. Because A is usually produced with a large velocity in collider experiments, the two decay photons will be 
largely boosted and seen by the detector (the electromagnetic calorimeter) as if it were a single-photon signal. 

To discuss the decay branching ratios of the other Higgs bosons, we need to specify all the parameters in the scalar 
sector of the THDM Lagrangian. In a CP-conserving THDM Lagrangian with natural flavor conservation (ensured 
by the discrete symmetry of <f>i — > ipi and </> 2 — > — ^2)7 there are eight parameters in its Higgs sector p5] |. They are 
mi, TO2, Ai,2,3,4,5, and ^12, or equivalently, Mh, Mh, M h +, Ma, a, tan/3, v, and /112. Of these eight free parameters, 
seven have been addressed in the previous section: four Higgs boson masses, two mixing angles (a and (3) and the 
vacuum expectation value v. There is yet another free parameter: the soft breaking ^12 term that so far has not been 
constrained.^ This is because up to the one loop order, ^12 does not contribute to the low energy data discussed 
above. With the assumption a = f3 — ir/2 we can write the hAA, HAA, HAZ and H + AW couplings as follows: 

C = l -\ h hAA + ^X H HAA + (24) 
Xhaz Z^Hd^A - Ad^H) + \ H+AW [W~(H + d^A - Ad»H+) + h.c] , 
where the coupling constants are given by 



6 The parameter /ii2 is defined through the interaction term — [if 2 'R'{4>i4>2}, which softly breaks the discrete symmetry of (j>i 
and </>2- It is given by 2A5U1U2 cos £ in Ref. it is sometimes written as m,3, see Ref. j3o| for example. 
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with M 2 = /Ltf 2 /(sin/3 cos/3). For a very light A, both the total width and the decay branching ratios of the other 
Higgs bosons in the model can be strongly modified to differ from the usual predictions of the THDM (with Ma at 
the weak scale v). 

At tree level, the partial decay width of h — > AA (or H — > AA) is given by 

A 2 

r(H)h->AA = ,n M • (25) 

It turns out that ft — ► is always the dominant (more than 90%) decay channel of ft except when 2M 2 ~ M 2 , in 
which case the hAA coupling is diminishing. When the parameter /112 increases, the decay width of the light Higgs 
boson can become large. For instance, for M% = 110 GeV and yUi2 = 200 GeV, Th—>AA — 34 GeV and 54 GeV for 
tan/3 = 1 and 2, respectively. Hence, in order for the considered model to be a valid effective theory we shall restrict 
the range of the parameter //12 so that the decay width of any Higgs boson should not be as large as its mass. For 
that reason, in the rest of this study, we shall consider the range of /112 to be < fi\2 < 200 GeV.[] In Fig. || we show 
the total decay width (r^) of ft, as a function of fix2 for a few values of Mh- Here, we take the value of tan (3 to be 0.5. 
The similar figure for tan/3 = 2 is identical to that for tan/3 = 0.5. This is because we are considering a — (3 — ir/2 
and the partial decay width of ft — > AA is unchanged after replacing tan /3 by cot f3. 

For Mh < 2mw, where raw is the mass of the gauge boson, ft — > bb is the sub-leading decay mode except when 
M 2 is in the vicinity of M 2 /2. For Mh > 2mw, the other decay modes (e.g., ft — > W + W~ , ZZ,tt) can open, and in 
that case, the AA mode is usually not the dominant decay mode. Since we are interested in the A — > 77 signal, we 
shall restrict our attention to values of Mh for which the decay branching ratio of ft — > AA can be sizable. To give a 
few examples, we show in Fig. || the branching ratios of ft — * AA , bb for Mh — 100, 115, 130 GeV with tan/3 = 0.5. 
(Again, the similar figure for tan/3 = 2 is identical to that for tan/3 = 0.5.) 

In addition to the AA channel, a heavy CP-even Higgs boson H can also decay into the AZ mode with a sizable 
branching ratio. The partial decay width for H — ► AZ is 

TH - AZ - ' (26) 



with Rzh — nv^/Mjj. (We note that Th^az does not depend on /ii2-) In Fig. 10, we show the total decay width of 
H as a function of /Ui2 for a few values of Mjj. Its decay branching ratios into the AA and AZ modes are shown in 
Fig. [ll] for various Mh values with tan/3 = 2. To study the tan/3 dependence, we also show in Fig. |l2[the branching 
ratios of H — > AA , ZA for /112 = 100 GeV. In this case, the sub-leading decay modes are H — > 66 , tt, etc. For M# 
less than twice of the top quark mass, the curves in Fig. |l^ are almost symmetric with respect to tan/3 = 1. This is 
again because we have set a — f3 — ir/2 and tan /3 is of order 1. 

The decay branching ratios oiH + are also largely altered in the THDM with a very light A, because the H + — > AW + 
channel becomes available. At the Born level, the partial decay width Th+->aw+ can be calculated using the formula 
given in Eq. (26) after substituting H — > H + and Z — » IV + Jj In Fig. [l3|, we show the dominant branching ratios of 
H + for the type-II model, in which the b — > S7 data requires M#+ > 165 GeV. 



7 For simplicity, we assume H12 to be a positive value, though the same conclusion also hold for a negative value of /ii2. 
8 The one loop corrections to T H +^ AW + have been calculated in Ref. they can modify the tree level width up to a few 
percent for the very low values of Ma considered here. 
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IV. PROBING A LIGHT A AT HIGH ENERGY COLLIDERS 



In this section, we discuss the potential of the present and future high energy colliders for detecting a light CP-odd 
scalar A to test our scenario of the THDM. 

The exciting feature of a light CP-odd scalar A is that the light CP-even Higgs boson h in the THDM can have a 
large decay branching ratio into the AA mode, and each A subsequently decays into a photon pair. For h with mass 
around 100 GeV, the decay particle A with mass around 0.2 GeV will be significantly boosted, so that the two decay 
photons from A are produced almost collinearly in the detector. When these two almost collinear photons cannot be 
resolved in the electromagnetic calorimeter, they will be reconstructed as a single photon. (The angular resolution for 
discriminating two photons in a typical detector will require Mu/Ma < 40.) As a result, the decay process h — ► AA 
will appear in the detector as a diphoton signature, and Z — > AAA as a triphoton signature, etc. Similarly, the final 
state of the production process e~e + — * ZAA will appear as a Z + 27 signature. 

In the following, we shall discuss in details the prediction of our scenario of the THDM on the decay branching 
ratio of Z — > AAA and the production rate of ZAA at the CERN LEP, Fermilab Tevatron, CERN Large Hadron 
Collider (LHC) and the future Linear Collider (LC). Other relevant production processes at the Tevatron and the 
LHC will also be discussed. Without losing generality, we again assume a — (3 — tt/2, motivated by the p-parameter 
constraint. Furthermore, we again take Ma to be 0.2 GeV such that the decay branching ration of A — * 77 is large 
and the life-time of A is short enough to be detected inside detectors of high energy collider experiments. 

A. The Decay Branching Ratio of Z —> AAA 

In the THDM, the Z boson can decay into the AAA mode via the two tree-level Feynman diagrams shown in 
Fig. [l4|. Since in the case of sm(/3 — a) = 1, the coupling of Z-A-h vanishes, only the diagram with the coupling of 
Z-A-H survives. For Mh > mz, due to the suppression factor from the three-body phase space, the partial decay 
width of Z — ► AAA, denoted as T(Z — > AAA), at the tree level, is small. For example, for Mh = 180 GeV, ^12 = 100 
GeV, and tan/3 = 0.5 or 2, T(Z — > AAA) = 1.32 x 10~ 8 GeV, which implies the decay branching ratio of the Z 
boson into three isolated photons (as identified by detectors) is of the order of 10 -8 and 10 -11 for tan/3 = 0.5 and 
2, respectively. (When Ma = 0.2 GeV the decay branching ratio Br(A — > 77) is 0.87 and 0.2 for tan/3 = 0.5 and 
2, respectively.) For a much heavier H, this tree level decay rate becomes negligible. In that case, a loop induced 
decay process might be more important. In Appendix |b[ we show the fermionic loop corrections to the decay width 
of Z — ► AAA, assuming that the other heavy Higgs bosons are so heavy that they decouple from the low energy data. 
We find that in general, this decay width is small unless the value of tan/? is very large. 

B. The Production Rate of ZAA at LEP and LC 

A light A could have been produced copiously at LEP-1 and LEP-2 experiments via the Yukawa process e~e + — > 
bbA Q. By searching for a light CP-odd Higgs boson in the associate production of the bottom quark pair, LEP 
experiments were able to exclude a range of Ma as a function of tan (3, when A decays into a fermion (lepton or quark) 
pair. On the contrary, the decay mode we are considering in this paper is A — > 77 which will likely register into 
detectors as a single-photon signal. Given the information on the decay branching ratio Br(A — > 77), it is possible to 
further constrain this model by examining the 667 events. However, as discussed in the previous sections, the (g — 2) M 
data has already constrained tan/3 to be small (less than about 2.6), so the production rate of e~e + — > bbA is not 
expected to be large at LEP. Here, we would like to consider another possible signal of a light A at LEP experiments, 
i.e. via the production process e~e + — > ZAA. 

The Feynman diagrams that contribute to the scattering process e~e + — > ZAA at the Born level are shown in 
Fig. |lj. With a — (3 — 7r/2, the tree level couplings Z-Z-H and Z-A-h vanish. Since by its definition, the mass of H 
is larger than that of h, the production cross section is dominated by the diagram (a) with h produced at resonance 
when Mh < — mz, where \/~S is the center-of-mass energy of the e~e + collider. Though the above observation is 
generally true, it is possible to have the value of //12 such that the coupling of h-A-A (i.e. Xh in Eq. ( p4| ) ) becomes so 
small that the production rate is instead dominated by the diagram with a H boson resonance. In that case, the event 
signature is to have a resonance structure in the invariant mass distribution of the Z boson and one of the A particles 
(i.e. one of the two photons observed by the detector), provided Mh < wS — Ma- Obviously, when Mh > yS— Ma, 
we do not expect any enhancement from the resonance structure, and the cross section becomes small. However, for 
a large value of [i\2, the width of h can become so large (cf. Fig. ||) that even for M h > — mz, the production 
rate of ZAA can still be sizable. The same effect also holds when the width of H becomes large. In Figs. |l6| and [Tt], 
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we show the production cross section for e~e + — > ZAA as a function of at the LEP and the LC for a few values 
of Mh and Mjj with tan/3 = 2 and tan/? = 0.5, respectively. (For completeness, we also give its squared amplitude 
in Appendix |(].) 

It is interesting to note that generally the whole complete set of diagrams for the scattering amplitude e~e + — * ZAA 
should be included in a calculation. For example, as shown in Fig. ^, when [iyi is about 58 GeV, the coupling 
vanishes for Mh = 130 GeV, and the bulk of the cross section comes from the diagrams (b) and (c) of Fig. [if]. 
Furthermore, the effect of interference among the complete set of diagrams can be so large that the distribution of 
the production cross section as a function of /X12 does not have a similar dip located the at the same value of [in . For 
example, as shown in Fig. |l7|, a broad dip in the distribution of <r{e + e'~ — > ZAA), predicted for LEP-2 with Mh = 130 
GeV (the dashed curve), is located at ~ 42 GeV, not 58 GeV. 

As noted previously, the experimental signature of the e~e + — * ZAA event is the associated production of a Z 
boson with two energetic photons. Based on this class of data sample, the LEP-2 "fcrmiophobic Higgs" search has 
imposed an upper limit on the decay branching ratio Br(/i — * 77) for a given Mh p2] ]. In that analysis, it was assumed 
that the production rate of e + e~ — * Zh is the same as the SM and the decay width of h is identical to that of the 
SM. One can express the experimental result on the photonic Higgs search in terms of the upper limit on the product 
a(e + e~ — > Zh) x Br(/i — > AA — > "77") as a function of Mh,. This upper limit can then constrain the allowed range of 
the parameter /U12 for a given Mh, which is shown in Fig. [TIL It turns out that the LEP-2 "fermiophobic Higgs" search 
data is not useful for constraining this model with tan/3 = 2. This is because in that case the decay branching ratio 
of A — ► 77 is about 0.1 (cf. Fig. 0) for Ma = 0.2 GeV, which largely reduces the Z + "77" event rate. Nevertheless, 
for tan/3 = 0.5, only a small region of /Z12 is allowed for Mh less than about 103 GeV.[] 

There is another important piece of data from LEP-2, that is the LEP SM Higgs boson search based on e + e~ — > 
Z — > Zh with h decaying into a 66 pair. It was reported |?3| that a handful events have been found to be compatible 
with the SM Higgs cross section for Mh about 115 GeV. Can a THDM with a light A be compatible with such an 
interpretation of the data? One trivial answer is to have Mh = 115 GeV with a choice of the free parameter /ii2 such 
that the decay branching ratio Br(/i — > 66) is about the same as that in the SM. This would obviously require the 
range of \in to be near the dip in Fig. ||. For a = (3 — n/2, /112 is about 53 GeV for a 115 GeV light CP-even Higgs 
boson h. This result implies a very specific production rate of the diphoton pair produced via gg — ► h — > AA at the 
Tevatron and the LHC. We shall come back to this production process in the following sections. 

Another solution to this question is to realize that the observed e + e~ — > Zh(—> 66) event rate at LEP-2 is determined 
by the product cr(e + e~ — > Zh) x Br(h — > 66). It can be the case that Mh is less than 115 GeV so that a(e + e~ — > Zh) is 
larger than that for the 115 GeV case. Because in a light A THDM the additional decay channel h — > AA is available, 
the decay branching ratio Br(/i — > 66) decreases. This reduction can compensate the increase in the production rate 
of Zh to describe the same experimental data. 

In Fig. [l8[ we show the corresponding range of the parameters /ii2 and Mh, assuming that the kinematic acceptances 
of the signal and the background events do not change largely as Mh varies.^] Hence, if we follow the LEP-2 conclusion 
on the SM Higgs boson search, i.e. at the 95% CL the current lower bound on Mh is about 113.5 GeV, then the 
product of cr(e + e" — ► Zh) x Br(h — > 66) cannot be larger than the SM prediction for Mh — 113.5 GeV. For a given 
Mh, Br(/i — > 66) cannot be too large. Therefore, this data could exclude the values of /ii2 near the dips shown in 
Fig. |^. As expected, this set of data and that for photonic Higgs search provide a complementary information on 
constraining the model. Combining these two sets of data, a light h with mass less than 103 GeV in the type-II 
THDM is excluded when tan/3 = 0.5. For Mh > 103 GeV, some constraints on the range of H12 can be obtained. The 
combined constraint is shown in Fig. [l^ for tan (5 = 0.5. For tan /3 = 1, a similar constraint can be obtained, and the 
region with Mh < 95 GeV is excluded. 

In the following analysis, we shall focus on the region of the parameter space that is consistent with the LEP Higgs 
boson search result, i.e. at the 95% CL the current lower mass bound on a SM Higgs boson is about 113.5 GeV, as 
discussed above. 



9 Based on the data presented in Ref. , this is the highest value of Mh for which we can obtain a useful bound on the value 

Of (1112. 

10 While this assumption is valid when Mh is around 115 GeV, it is likely to fail when Mh is close to the Z-boson mass (mz)- 
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C. The Production Rate of ZAA and W AA at Tevatron and LHC 



At the Run-2 of the Fermilab Tevatron, a 2 TeV proton-antiproton collider, the production rate of pp — > ZAA is 
about 0.1 pb for M h = 110 GeV, /Lt 12 = 100 GeV, M H = 1 TeV, and tan/3 = 0.5 or 2. Being a hadron collider, 
Tevatron can also produce a light A pair via the constituent process q'q — > W ± AA. The Feynman diagrams for this 
scattering process are the same as those depicted in Fig. [15] after replacing Z by VF ± everywhere and H by ff* in 
diagram (b). For the same parameters given above, the W^AA production rates is about 0.2 pb. For completeness, 
we show in Figs. [l9] and ^0| the cross section for ZAA and W^AA productions, respectively, at the Run-2 of the 
Tevatron, and the LHC (a 14 TeV proton- proton collider). It is interesting to note that when the mass of the charged 
Higgs boson is not too large (consistent with the case of tan (3 = 2), and the coupling of h-A-A vanishes, the production 
rate of W^AA is dominated by the associate production of A and i? ± which subsequently decays into a W^-boson 
and A. Since the experimental signal of A is an "isolated photon" , this signal event appears as an event with a 
and two photons, hence, its SM background rate is expected to be small. For completeness, the squared amplitude 
for this partonic process is also presented in Appendix ^|. 

D. Other Production Mechanisms of a Light A at Colliders 

In addition to the above production processes, a light A can also be copiously produced at hadron colliders, such 
as Tevatron and LHC, via qq,gg — > bbAorttA, gg_—>h—> AA, gg — > H — > AA, and qq — > H+H~ — > W + W~AA. 
Because of the potentially large background, the bbA mode is not likely to be observable. However, the ttA mode can 
be easily identified by requiring an isolated photon with a large transverse momentum. At the Tevatron Run-2, the 
inclusive rate of ttA with a 175 GeV top quark is 9.8(0.6) fb for tan/3 = 0.5(2), and at the LHC, it is 1.6(0.1) pb. The 
production rates for the last two processes can be easily calculated by multiplying the known cross sections for the 
production of gg — > hor H [pMj and qq — > H + H~ p5| by the relevant decay branching ratios (given in the previous 
sections). For example, at the Tevatron Run-2, the production cross section of gg — ► h — > AA is 0.8 and 0.5 pb for 
Mh = 110 GeV and 130 GeV, respectively, with /L«i2 = 100 GeV and tan/3 = 0.5 or 2. Since we have set a = j3 — ir/2 
(motivated by the /^-parameter constraint) in all our calculations, the production rate of gg — ► h is independent of 
tan/3. Furthermore, when //12 = 100 GeV, the decay branching ratio for h — > AA is about 1, cf. Fig. ||. Therefore, 
the above rates for tan/3 = 0.5 or 2 are about the same. At the LHC, the rates are 35 pb and 27 pb for Mh = 110 
GeV and 130 GeV, respectively. Hence, this mode is important for further testing the THDM with a very light A. In 
Fig. pTL we show the production rate of gg — ► h — > AA at the Tevatron Run-2 and the LHC as a function of /i 12 for 
various 

The signal rate of gg — ► H — > AA is different from that of gg — ► h — > AA because the decay branching ratio of 
H — > AA is not the same as that of h — ► AA and the coupling of H to £ in the loop has a factor of cot (3, cf. Eq. Jlc|). 
For instance, the branching ratio of -> AA is about 0.12 (0.47) for M H = 180 GeV (350 GeV) with fi 12 = 100 GeV, 
so that the cross section of gg — > H — > AA is about 6 fb (2.6 fb) at the Tevatron, and 0.48 pb (0.89 pb) at the LHC 
when tan/3 = 2. For tan/3 = 0.5, the low energy data requires Mh to be around 1 TeV, so that its rate is negligible. 
In Fig. |22j, we show the pp,pp(gg) — > H — > ylvl production cross section as a function of ^12 at the Tevatron Run-2 
and the LHC for a few values of Mh, with Ma = 0.2 GeV and tan/3 = 2. Note that when the rate of gg — > i? — * 
is small (for certain values of /X12), the rate of gg —> H —> AZ becomes large because the sum of the decay branching 



ratios of these two modes is about 1, cf. Fig. 11, Hence, their roles to the discovery of a heavy Higgs boson in this 
model are complementary to each other. 

Usually, a charged Higgs boson H + in the THDM is assumed to decay via the heavy fermion pairs, either the tb, 
cb or t + v modes. However, as shown in Fig. |l3|, when A is light, the decay mode of W + A can become dominant. In 
that case, the scattering process qq — * H + H~ — > W + W~ AA will be seen by the detector as a VF-boson pair with 
two isolated photons. With a proper kinematic cut, this event can be separated from its SM backgrounds. (Although 
we have limited ourselves to the discussion on the decay mode of A — > 77, the other decay mode into an e + e~ pair 
can also prove to be useful for testing such a model.) For M H + = 180 GeV (350 GeV), the decay branching ratio of 
H~ -> W~A is about 1 (0.8), so that the cross section of qq -> H+R- -> W + W~AA is about 1.7 fb (0.01 fb) at 
the Tevatron, and 31 fb (1.5 fb) at the LHC when tan/3 = 2. Again, when tan/3 = 0.5, the low energy data requires 
M H + to be around 1 TeV, so that its rate is negligible. 



11 To compare with the experimental data, one should also include the decay branching ratio Br(A — > 77), cf. Fig. [?], for each 
CP-odd scalar A decaying into its photon mode. 
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In the THDM, there is no tree level coupling Z-W^-H^ , therefore the cross section for the scattering process 
qq — > W ± H T — > W + W~A at the tree level is dominated by the bottom quark fusion bb — > W ± H T — > W + W~A. At 
the Tevatron, its rate is negligible (about 0.1 fb), and at the LHC, its rate is about 0.1 pb for M H + = 200 GeV, and 
tan/3 = 2. 

At the future Linear Collider (a 500 GeV e + e~ collider), the ZAA production rate is shown in Figs. |l^ and 
[Tt| for a few choices of parameters. Besides this production mode, a light A can also be copiously produced via 
e+e~ — > H + H~ — > W + W~~ AA provided that the cross section e+e~ — > H + H~ is not small. For Af H + = 180 GeV 
(200 GeV), we expect that rate of W+W'AA to be about 41 fb (22 fb) for tan/3 = 2. (The Br(iJ+ -> yW+) is 
about 1 and 0.9 for M H + = 180 GeV and 200 GeV, respectively.) 

V. DISCUSSION AND CONCLUSION 



In this paper, we examined the possibility of having a very light CP-odd scalar A (with a mass about 0.2 GeV) in 
a general THDM. After examining the relevant low energy data, we found that this model is either excluded already 
(according to the DH prediction of the SM contribution to the muon (g — 2) M ) or its parameter space has been largely 
constrained. Assuming the CLY prediction of the SM contribution, the muon (g — 2)p data requires tan/3 < 2.6, 
regardless of the other parameters of the type-II THDM. (For the type-I THDM, the muon (g — 2)„ data requires 
tan/3 > 0.4.) For such a light A, the CUSB data on Br(T — > Aj) requires tan/3 to be around 1. For a type-II THDM, 
the b — > S7 data requires M H + to be larger than about 165 GeV to 200 GeV when tan/3 is larger than 1. For a 
type-I THDM, the constraint on M H + is much looser. The p-parameter data also imposed a stringent constraint on 
the difference between Mh and M H + for a given value of sin(/3 — a). For sin(/3 — a) — 1, Mh and M H + have to be 
almost equal. The Rf, data also provides a stringent bound on the model. For example, when tan/3 = 0.5, the i?f, data 
requires M H + to be around 1 TeV in the THDM. Consequently, due to the p-parameter constraint, Mh should also 
be around 1 TeV. A summary of the low energy constraints on this model is given in Figs, [l] - || as well as Tabic B. 

After finding the allowed parameter space of the model, we examine the impact on the decay branching ratios 
and total decay width of the Higgs bosons due to the presence of a light A. Depending on the value of the soft- 
breaking parameter /112, present in the Higgs potential of a general THDM, the total decay width of h, H or H + 
can become large because of the large phase space volume for the decay channels (h — > AA,AZ), (H — > AA,AZ), 
or (H + — > AVF + ). To have a valid perturbative calculation, we only consider the values of ^12 such that the total 
decay width of the Higgs bosons is small as compared to its mass. Due to the small mass of A, the decay branching 
ratios for the decay modes AA, AZ, or AW^ can be sizable, which can result in a very different detection mode for 
the THDM. In Figs. §-0, we showed a few of such examples. 

The exciting feature of such a light A is that when it decays into a photon pair, because of the typical large energy 
of A produced from the decay of other heavy Higgs bosons, its decay photon pair will register in the detectors as a 
single photon signature. Hence, the SM background rate for detecting such a signal event is expected to be generally 
small. Therefore, the Tevatron Run-2, the LHC and the future LC have a great potential to either detect a light A 
in the THDM or to exclude such a theory model. A few potential discovery modes at various colliders were given in 
Section IV, cf. Figs. |l§|2|. 
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Note Added: 

During the preparation of this manuscript, two similar papers |39[| were posted to the hep-ph archive very recently. 
For the part we overlap, our results agree. Furthermore, after the submission of this paper, a new article pCj ] 
concluded a tighter bound on M H + from b — > sj data than the one [M adopted in our analysis. Nevertheless, the 
general conclusion about the phenomenology of a light A discussed in this paper remains unchanged. 
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APPENDIX A: ANOMALOUS MAGNETIC MOMENT OF MUON IN THE THDM 



Given the THDM interaction Lagrangian 



£ CsP-^S + iC A flj 5 ^A+ (y2C H ±lLv L B- + h.c.H , (A.l) 

\S=H,h J 



2lTL\y 

the contribution from the diagrams depicted in Figs. |23| (a) and (b) reads 



, a m 2 

one — loop H_ 



- X / >I>-C: <n-r). (A.2) 

where j is summed over h, H, A, and H . The respective g\x) are: 

x 2 (2 - x) 

c 2 +h H A l - x Y 



9 H ' h (*) = 5 X \ (2 t , (A.3) 



^ = ^TMM ' (A - 4) 
^<*> = iffe^T (A - 5) 

with Aj = (Mj/m^) 2 . The respective Cj are given in Table g for a type-I or type-II THDM. 

As for the contribution from the Barr-Zee diagram depicted in Fig. 2^ (c), the dominant two- loop contribution 
comes from a CP-odd scalar in the loop. There is another diagram with the virtual photon replaced by a Z gauge 
boson. Because it is highly suppressed, we will not take it into account. According to [fl5|| , the dominant two-loop 
contribution to a p is given by 



a two- 10 op = « 2 rn£t £ N{Q}C fXf f( Xf ), (A.6) 

87T S w M w f=t br 

with Xf — ( m f/MA) 2 , N£ — 1 (3) for leptons (quarks), to/ and Qf are the mass and charge of fermion, and Cf is 
given by the interaction Lagrangian of the CP-odd scalar to fermions: 

C = i 9 -^h 5 fA, (A.7) 

where C t = cot/3 for either a type-I or type-II THDM, whereas Cf — — cot/3 and tan/3 for a type-I and type-II 
THDM, respectively, for / = b, n, and r. Finally, f{x) is given by 



I) 



x — z(l — x) 



For comparison, we note that if only one-loop contributions are included in our analysis, we will derive a different 
constraint on the parameters of tan/3 and Ma- In that case, the figures corresponding to Figs. || and § are shown as 
Figs. p3and|25|. 



APPENDIX B: THE FERMION-LOOP CONTRIBUTION TO THE Z -> AAA DECAY 

Under the scenario that the branching ratio of A — > 77 is close to unity, it is possible that the CP-odd scalar A 
has a large coupling to fermions only. This can happen when the other scalars do not exist at all or are so heavy 
that they would decouple from the low energy effective theory because of their small couplings to A. We consider 
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this assumption to examine the contribution of fermion loops to the decay Z — » A^4A[^] One of the fermion loop 
diagrams contributing to the Z(p) — > A(ki)A(k2)A(k3) process is shown in Fig. |26|. (There are five other diagrams 
with the obvious permutations of the pseudoscalar boson momenta.) 

The fermion loop contribution to the decay of Z — > AAA was first roughly estimated in |[(|, and then was re- 
examined in fl37[ by considering only the top quark loop contribution. A complete calculation, including also bottom 
quark contribution with a large bbA coupling, was never presented in the literature. We consider an effective theory 
in which the coupling of the CP-odd scalar to quarks is C q m q /v, where v(= 2mw/g) is the vacuum expectation 
value, m q is the mass of the quark, and the coefficient C q depends on the choice of models. For a SM-like coupling, 
assuming the existence of A, C q — 1 for both up- and down-type quarks. For a type-I or type-11 THDM, C q — cot (3 
for up-type quarks. For down-type quarks, C q — — cot j3 and tan (3 for a type-I and type-II THDM, respectively, 
Since the coupling of Z to up- and down-type quarks is fixed by the SU(2) L x U(l)y gauge interaction, the effective 
Lagrangian can be written as 



C = J2 



j_ 

2c, 



T, <v a (n q rfl ^ 5 )n 7 , id m qCq - „,5 

~ u ql \9v ~ 9a1 ) u q H -z 

, Zmw 



U q 7°U q A 



(B.l) 



Given the above interaction, the Z — > AAA amplitude can be expressed, in terms of Mandelstam variables s 
(ki + fc 2 ) 2 , t = (fc 2 + fc 3 ) 2 , and u = (fci + k 3 ) 2 , as 



M = 



where 



F{t,u)k? + F{u,t)k% e Q (p), 



(B.2) 



F(t, u)=Y, C «[^r) [^ M a ~ 3 - *)(fo(a, u) - D (t, u)) + (t - s)D (s, t) 

8 \ m Z J 

+ 4(C (t) - CoOO) 



(B.3) 



with the scalar three- and four-point integrals given by 



l-K 1 J 



[k] [k-k 3 ] [k-p]' 



(B.4) 



D (s,t) 



d 4 k 



in 2 J [k] [k - fci] [fc - fci - fc 2 ] [fc - p] ' 



(B.5) 



where we use the short-hand notation [£] 

the pseudoscalar bosons momenta. The squared amplitude, after averaging over the Z boson polarizations, is 



i 2 +ie. The remaining scalar integrals can be obtained by permuting 



{{m z - M A ) 2 - t){{m z + M A ) 2 - t)\F{t, u)\ 2 



+ ({m 2 z - M 2 A ){t + u) + M% + tu) Re(F(t, u)F\u,t)) + (t u) 



(B.6) 



Given this result, the decay width can be computed by the usual methods. Although our result is quite general, we 
only expect important contributions from both b and t quark loops: contributions from lighter quarks are suppressed 
by the factor {m q /mz) A and can be ignored unless C q is extremely large, which however is unlikely to be true because 
of the tree-level unitarity constraint on the couplings. As we are primarily interested in examining the case of a very 
light pseudoscalar Higgs boson, we can safely neglect M A in our calculation. In the M A — > limit the scalar integrals 
become |3ql 



C (s) 




4m 2 



(B.7) 



The decay mode Z — > AA is forbidden by the Yang-Landau theorem. 
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D (s, t) = - [G(a, t, s) + G(s, t, t) - G(s, t, m|)] , 
where F is given in Eq. (p3|), and G can be written in terms of Spence functions 



(B. 



G(s,t,y) = - i Sp 



Sp 



- b 



-Sp 



a + — b 



Sp 



log ( — ) log ( 1 - it 



b — a+ 

y 



(B.9) 
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with (A = . / 1 



4m2 (m| - s-i)/(s t), a± = |(l±</»)/2, and 6 = ±(l + y^l - 4(m% - ie)/y) . 

It is easy to estimate the order of the scalar integrals arising from a heavy quark: if m q > 2mz and Ma — > 
0, we can use the heavy mass expansion approximation. The leading term of the three-point scalar integral is 
C(s,t) ~ — l/(2m 2 ), which can be differentiated with respect to m 2 to give D(s,t) ~ l/(6rn 4 ). We then have 

C(s) 1.62 x 1(T 5 GeV~ 2 and D(s,t) ~ 1.78 x 1CT 10 GeV~ 4 for m q = 175 GeV, the top quark case. On the 

other hand, for the b quark loop, numerical evaluation gives | C(a) | ~ 10~ 4 — 10~ 3 GcV~ 2 and \D(s, t)\ ~ 10~ 6 — 10~ 5 
GeV -4 , which indicates that bottom quark contributions may compete with that of the top quark, even if we consider 
the factor of {m q /mz) 4 in (B.3). In conclusion, the decay branching ratio due to the top and bottom quark loop 
contributions is 



Br(Z -> AAA) = 1.3 x KT 18 C? + 2.47 x l(T 17 C b 6 + 7.63 x l(T 18 C t 3 C & 3 



(B.10) 



which is many orders below the estimate given before. However, there is no contradiction because the authors of |37|] 
used a rough estimate to show that top quark loops cannot enhance the branching fraction of Z — > AAA beyond 10^^ 
in the case of Ct — 1. From (B.10) we can see that top quark contribution is smaller unless C* 3> 1, while bottom 
quark contribution is larger for C& > Ct- A large Cb in a type-II THDM-like model implies large tan/3; however, 
because of the unitarity bound, the coupling Cb cannot be arbitrary large. By requiring the validity of a perturbation 
calculation, we can derive the upper bound on Cb to be about 120 which yields Br(Z — > AAA) ~ 10~ 5 . 



APPENDIX C: THE PROCESSES // -> ZAA AND qq' -> W + AA 
1. Squared Amplitude of // — > ZAA 



The Feynman diagrams contributing to the scattering of // — ► ZAA are shown in Fig. |T^. The scattering amplitude 
for /(pi)/(pa) — * Z{k)A{k\)A{k<x) can be written as: 

M = u flu Up l + g f R P R ) u } \F Q + F 1 k$k% + F 2 k£k%\ e* v (k), 



(C.l) 

where g^ L = 1$ — ejs^ and g R = —efS^, with I3 the weak isospin of the fermion and ej its electric charge in units 
of positron's. After averaging over the spins and colors of the initial state and summing over the polarizations of the 
final state particles, we obtain 



iEi^i 2 = 



N c g 6 X 2 



spins 



4m| (s — m 2 z ) 2 



2 2 

*ij=0, j>i 



i=0 



(C.2) 



where N c — 1 or 1/3 for leptons or quarks, respectively. Furthermore 

2 if - 2e/J 3 «* + 2e 2 s 4 

V = - 



(C.3) 



whereas the Q are 



Co = 4(fc -pi)(k -p 2 ) +m 2 z s, 



(C4a) 



Ci = {M\ s - 4(fc 2 - Pl )(k 2 .pa)) {M\m\ - (k ■ k x f) , 



(C4b) 
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C 2 = Ci(fci^fc 2 ), (C4c) 

Coi - 2 (m| v(h) ~ (fc • h) r,(k)) , (C4d) 

Co2-Coi(fci^fc 2 ), (C4e) 

Cia = 2 (m|(fci • k 2 ) - (k ■ fc x )(fc ■ fc 2 )) r^), (C4f) 

j 1 {i) = s{l-k 2 )-2({l-p l )(k 2 ■ p 2 ) + {I ■ p 2 )(k 2 - pi)), (C4g) 
and s = (pi +p 2 ) 2 . The form factors £j are given by 

& = = (c.5) 

Diagrams with the ZZAA and ZZH(h) vertices contribute to the £o form factor, whereas those with ZAH(h) vertex 
contributes to £1.2 



1 1 I ~***l> 



(2c s+ ^l 2 - s%pS S _M\ - s 2l3 (c a c 3 p - s a s 3 ^ Mlj (f) = h, 



g'- ow-'q U,-*, <C6b> 

where <5± = /3 ± a and S2 = {k\ + fc2) 2 . In addition, it is understood that the Higgs propagators acquire an imaginary 
part in the resonance region, i.e. p 2 — M? — > p 2 — A/ 2 + iF^Mj,, where is the total width of <f>. 

2. Squared Amplitude of qq — > W + AA 

The partonic process gg' — > W + AA receives contributions from just five diagrams if the quark masses are neglected. 
These diagrams can be obtained from those contributing to the process // — > ZAA, cf. Fig. [l5|. After a few changes, 
t he a bove results can also be easily translated to obtain the respective squared amplitude. First of all, in equation 
( p.lD we have gh — Vqq' /2 and gn — 0, together with the approp riat e change of notation regarding the Dirac spinors. 



V qq i is the CKM mixing matrix element. Secondly, in equation (C.2) \ is now given by 



X 2 = (C7) 
whereas the substitution mz —> myy must be done everywhere. Finally, the £j form factors are defined now as 

& = = # +ef + ef + , (c.8) 



where £g are the same as those in ( C6a ), and 

^ = (HM 2 -4 + ' (r ' ! " 



Again, in the resonance region, the charged Higgs boson propagator acquires an imaginary part. 
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FIG. 1. The regions (below the curves for type-II and above the curves for type-I THDM) in the tan/3 versus Ma plane 
allowed by the data at the 95% CL. Three different curves are displayed depending whether the SM prediction is given by 
the CLY, AY or J calculation. There is no allowed region in this range of parameters according to the DH calculation. (Here, 
a two-loop calculation for the THDM contribution, cf. Appendix [X], is used.) 
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FIG. 2. The regions (above the curves) in the tan /3 versus Ma plane of a type-I THDM allowed by the a M data at the 95% 
CL. There is no allowed region in this range of parameters according to the DH calculation. (Here, a two-loop calculation for 
the THDM contribution, cf. Appendix^, is used.) 
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FIG. 3. The regions in the tan/? versus Ma plane of a type-II THDM allowed by the a M data at the 95% CL. The allowed 
regions based on the calculations of CLY, AY and J are below the curves. The region allowed by the DH calculation is bounded 
by the dotted line. (Here, a two-loop calculation for the THDM contribution, cf. Appendix |a| is used.) 
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FIG. 5. Allowed parameter region for a THDM with a light pseudoscalar A. The upper bounds on tan (3 derived from the 
a M data is for a type-II THDM with M A = 0.2 GeV. The lower bound from R b holds for either type-I or type-II THDM (when 
tan/3 ~ 1). The b — > s-y data does not provide any useful constraint for a type-I model when tan/3 > 1. 
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FIG. 13. The H + -> AW + and H+ -> tb branching ratios for two values of tan (3 in the type-II THDM. 
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FIG. 16. ZAA production cross section at LEP-2 and LC with yS = 210 GeV and 500 GeV, respectively, for a few values 
of M h and M H with M A = 0.2 GeV and tan/3 = 2. 
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FIG. 17. ZAA production cross section at LEP-2 and LC with \/~S = 210 GeV and 500 GeV, respectively, for a few values 
of M h with M A = 0.2 GeV, M H = 1 TeV and tan/3 = 0.5. 
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FIG. 19. ZAA production cross section at the Run-2 of the Tevatron and the LHC, for a few values of Mh, with Mh = 1 
TeV, M A = 0.2 GeV and tan/3 = 0.5. 
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FIG. 20. W ± AA production cross section at the Run-2 of the Tevatron and the LHC, for a few values of Mu, with M H ± = 1 
TeV, M A = 0.2 GeV and tan/3 = 0.5. 
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FIG. 21. The pp,pp(gg) — > h — » AA production cross section as a function of /112 at the Tevatron Run-2 with y/S = 2 TeV 
and the LHC with y/S = 14 TeV, for a few values of M h with M A = 0.2 GeV and tan/3 = 0.5 or 2. 
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FIG. 22. The pp,pp(gg) — ► H — > production cross section as a function of /ii2 at the Tevatron Run-2 with y/S = 2 TeV 
and the LHC with \/S = 14 TeV, for a few values of M H , with M A = 0.2 GeV and tan/3 = 2. 
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FIG. 23. Contribution from the THDM to the anomalous magnetic moment of muon: neutral Higgs bosons (a), charged 
Higgs boson (b), and the leading two- loop contribution from the CP-odd scalar (c). 




FIG. 24. The regions (above the curves) in the tan (5 versus Ma plane of a type-I THDM allowed by the a M data at the 95% 
CL, based on a one-loop calculation. There is no allowed region in this range of parameters according to the DH calculation. 
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FIG. 25. The regions (below the curves) in the tan f3 versus Ma plane of a type-II THDM allowed by the a M data at the 95% 
CL, based on a one-loop calculation. There is no allowed region in this range of parameters according to the DH calculation. 




FIG. 26. Representative fermion loop diagram to the Z — > AAA decay. 



Constraint 


Type-I THDM 


Type-II THDM 


(s-2)„ 


tan/3 > 0.4 


tan/3 < 2.6 


(tan 13 > 1) b —> S~f 


M H+ > 100 GeV 


M H + > 200 GeV 


(0.5 < tan < 1) b — > S7 




M H+ > 200 - 350 GeV 


(0.6 < tan ^ < 1) Rf, 


M H + > 200 - 600 GeV 


M H + > 200 - 600 GeV 


( S raQ3 - a) = 1) Ap 


M H ~ M H + 


M H ~ M H + 


(sta(/3 - a) = 0.8) Ap 


M H ~ 1.2 M H + 


M H ~ 1.2 Mu+ 


(staOB - a) = 0.5) Ap 


M„ ~ 1.7 M H + 


M H ~ 1.7Af H+ 



TABLE I. Constraints from the low energy data for type-I and type-II models, with Ma = 0.2 GeV. Here, we use the CLY 
calculation for the SM prediction of a,j. When sin(/3 — a) is fixed to be 1, there is no dependence in p, otherwise, we assume 
M h = 110 GeV. 



E 1 1 1 1 1 = 

: THDM-ll : 
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THDM 


C H 


c h 


Ca 




I 


sin a/ sin {3 


cos a/ sin /3 


cot/3 


cot/3 


II 


cos a/ cos/3 


— sin a/ cos /3 


— tan/3 


— tan/3 



TABLE II. Higgs scalars-muon couplings in the THDM. 
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